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22 Abstract 

23 In this paper the realized niche of the Snowy Plover (Charadnus alexandrinus) , a 

24 primarily resident Florida shorebird, is described as a function of the scenopoetic and 

25 bionomic variables at the nest-, landscape-, and regional-scale. We identified some pos- 

26 sible geomorphological controls that influence nest-site selection and survival using data 

27 collected along the Florida Gulf coast. In particular we focused on the effects of beach 

28 replenishment interventions on the Snowy Plover (SP), and on the migratory Piping 

29 Plover (PP) {Charadnus melodus) and Red Knot (RK) (Calidris canutus). Addition- 

30 ally we investigated the potential differences between the SP breeding and wintering 

31 distributions using only regional-scale physiognomic variables and the recorded occur- 

32 rences. To quantify the relationship between past renourishment projects and shorebird 

33 species we used a Monte Carlo procedure to sample from the posterior distribution of the 

34 binomial probabilities that a region is not a nesting or a wintering ground conditional 

35 on the occurrence of a beach replenishment intervention in the same and the previous 

36 year. The results indicate that it was 2.3, 3.1, and 0.8 times more likely that a region 

37 was not a wintering ground following a year with a renourishment intervention for the 
3a SP, PP and RK respectively. For the SP it was 2.5. times more likely that a region 

39 was not a breeding ground after a renourishment event. Through a maximum entropy 

40 principle model we observed small differences in the habitat use of the SP during the 

41 breeding and the wintering season. However the habitats where RK was observed ap- 

42 peared quite different. While ecological niche models at the macro-scale are useful for 

43 determining habitat suitability ranges, the characterization of the species' local niche 

44 is fundamentally important for adopting concrete multispecies management scenarios. 

45 Maintaining and creating optimal suitable habitats for SP characterized by sparse low 

46 vegetation in the foredunes areas, and uneven/low-slope beach surfaces, is the proposed 

47 conservation scenario to convert anthropic beach restorations and SP populations into 
4a a positive feedback without impacting other threatened shorebird species. 



1 Introduction 



The increasing availability of spatio-temporal data on species presence, along with the avail- 
ability of remotely sensed data and GIS techniques, has greatly enhanced in the last decad e 



2006: 



Soberon 



2007) 



52 the study of the distribution of thousands of species (lElith et al 

53 However, the individuation of the species' range is performed in a Grinnellian way, con- 

54 sidering mostly scenopoetic variables that are suitable to describe the fundamental niche 



55 ( Soberon 



20071 ). Studies on the distribution of species require the consideration of biotic 



variables (Eltonian perspective) in order to truly characterize the realized and the funda- 
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mental niches (jColwell and Range! 120091) . If the habitat exhibits conditions that lie entirely 
within a species' niche, a population persists without immigration from the external world, 
whereas if conditions lie outside the niche, the species faces possible extinction. Analysis 
of species' niches are essential to understand controls on species' geographical range limits 



and how these limits might shift in response to climatic changes (Holt 



2009; 



Zimmermann et al 



2009; 



Tinglev et al 



20091) . The Hutchinson' s duality consists in considering simul- 
taneously exogenous variables that describe the biotope in which the species live, and the 
biotic variables that characterize the inte ractions of the species with other living and non- 



and landscape ecology 



living controls (jColwell and Range! _ 
( Knowlel 2OO9I 



2009) . Recentl y, the emerging hclds of phylogcography 



Wang . 



20101) have significantly improved species dis- 



tribution modeling and to detect d i fferences among species inc l uding data at the ce l l-level 



(e.g. DNA sequences) ( 


Fu 


Qk et al. 




200' 


7 : 


Rissler and Apodaca 


Kearnev and Porter 


20091 


Miller et al. 




2009|). 



2007 



K'upper et al 



2009; 



71 The selection of breeding and wintering habitat by shorcbirds and their consequent sur 

72 vival may be influenced by a combination of factors, includin g human recreational activities 

73 predator activ i ty, prey availability, and the habita t substrate (jHoover and Brittingham . 



Newton 



1998 



Jones 



2001 



Colwell et al 



1998 



2007a|) . Nest-site selection and nest-survival pat- 



terns reveal in general an influence by a combination of the aforementioned environmental 
and biological factors working in concert in addition to physical features surrounding the 
nest-site. Few avian habitat studies have been able to compare multipl e ecological h ypothc- 



ses of species distribution and multispecies nest-site selection decisions (jJones 



200 ll) to aide 



management policies. While habitat selection is ofte n assumed to be adaptiv e , evidence fo r 



1999; 



Jones 



2001) 



adaptive habitat selection in birds has been mixed (jClark and Shutler . 
The consideration of multiple predictors collectively for detecting species distribution, is use- 
ful for habitat management, and it benefits the conservation of rare and declining species. 
Shorcbirds reproductive success is correlated with the stability and quality of the nesting 
environment. In particular, here we study the effects of beach replenishment (or renour- 
ishmcnt) on the Snowy Plover (Charadrius alexandrinus) , a state-threatened shorebird in 
Florida (Figure [TJ. SP females in general show fidelity to nesting beaches, making artificial 
beach nourishment practices and the subsequent physical and biological changes to habitat 
directly relevant to their recovery. The reduction in reproductive output is in general pri- 
marily a consequence of decreased nesting success. The result of reduced nesting success 
is more precisely described as reduced juvenile recruitment rather than reduced number of 
nests, since nesting attempts still occur but are not successful. However late renourishments 
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92 in the wintering season can be a source of disturbance for the Snowy Plover. In comparison 

93 to the SP, two other shorcbird species have been considered. The Piping Plover, Charadrius 

94 melodus (federally designated as threatened), and the Red Knot, Calidris canutus (threat- 

95 ened in New Jersey, and a candidate for Endangered Species Act protection), are migratory 

96 shorebirds w hose wintering/stopover time in Florida is on av erage 3 months and 3 weeks 



97 respectively (jHarrington . 



2001 



Elliott Smith and Haig 



20041) (Figure [TJ. The Red Knot 



98 subspecies Calidris canutus rufa is the only endangered species among the species consid- 

99 ered, under the federal Endangered Species Act. It has been established that the Calidris 

100 canutus rufa uses some Florida beaches as stopover areas during its migratory route to South 

101 America. The Wilson's Plover (WP) is an other resident shorebird reputed to be the main 

102 competitor of the SP, however it is not considered in this study due to its least-concern status 

103 in Florida. 



Due to habitat loss, many threatened, endangered, and at-risk species (TER-S) are de- 
creasing in abundance. The a-priori evaluation of the effect of beach restoration activity 
on species is fundamentally important to understanding the effectiveness of the intervention 
and to optimizing strategics. Beach reno urishment is mainly carried out to preserve existing 



structures and to increase the beach area ( Smith et al 



20091 ) . This potentially translates into 



new income from tourism and beac h activ i ties. I n Florida, an average of $ 90 million is spent 



annually on beach renourishment (jPSDS 



20101 ). Over the past two decades, more than 50 



n2 large renourishment projects have been undertaken in the state, with a typical project aver- 
n3 aging approximately 4.5 km in length. More than 242 km of G ulf and Atlantic coa st beaches 



have been impacted by renourishment sand during that time (jWang et al 



20051) . The five 



Gulf st ates account for more than forty percent of all renouris hment 



States ( PSDS 



activi ty in the United 



19961 ). Beach renour- 



20101) . Florida alone accounts for thirty percent flFinkll . 

n7 ishment efforts are not without significant social or legal impacts. In June 2010 the Supreme 

us Court handed down a unanimous verdict that effectively allows the Florida state government 

n9 to resume beach renourishment projects without paying for property that homeowners claims 

120 have been "taken" . In 2003 a group of NW Florida coastal homeowners protested against a 

121 replenishment, claiming lowered property values due to the increase in beach width (on aver- 

122 age 50 m) and subsequent greater public access to the beach. The US Supreme Court rejected 

123 the appeal, declaring the renourishment as a necessary intervention for preserving the coastal 

124 ecological communities and human structures especially in ligh t of the increase in sea-level rise 



125 and of extreme meteorological events due to climate change (ISC-USA 



2010T ) . As a societal 



issue, beach renourishment is one of the most expensive interventions in civil and environmen- 
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tal engineering, considering also the environmental variables that are often unpredictable and 
the factors affecting the coastline due to c limate change and e xtreme climatological events 



(Smith et al 



2009; 



Gopalakrishnan et al 



2010; 



Landry 



20101 ) . For instance unpredictable 



variations in ocean energy impact due to littoral currents, or strong hurricanes, can rapidly 
destroy the renourished areas. Renourishment pr ojects are rare ly designed to incorporate the 



2010b() . The primary considerations 



life-history needs of shoreline-dependent species (BBCS, 
in planning a renourishment are sand source, and compatibility of the borrowed sand with 
the native beach, including grain size, composition and color. Disturbances associated with 
beach replenishment, such as dredging and sand-pipeline movement, represent an additional 
and potentially significant barrier to breeding and nesting for both shorebirds and water- 
birds. For example, investigators have recommended avoiding beach management practices 
that disturb beach microhabitats (e.g., epheme ral pools and bay tidal flats) importa nt for 



2000; 



Grippo et al 



2007 ). Re- 



Snowy Plover and Piping Plover chick survival (jElias et al. 
plenishment entails substantial changes in beach morphology that potentially cause changes 
in local movement patterns, resting behavior, or habitat use of shorebird species during 
the tidal cycle. Similarly, the potential disturbance to benthic macroinvertebrate assem- 
blages could alter fee d ing behavior in bird specie s whose diet relies on benth ic organisms 



([Bishop and Peterson . 



2005 



Duean and Hubbarc 



2006 



Peterson et al 



2006). Beach re- 



us nourishment has been found to alter shorebird distributions more than scabird distributions 



(jGrippo et al 



20071 ). However not many studies exist o n the topic and the effect of re 



2007) 



Lottl (120091 ) 



nourishments seems to be very strongly species-dependent (jGrippo et al 
analyzed the effect of sand replacement projects on SP and PP along the F lorida beaches 
report ing a qualitative negative correlation, without indication of the causes 



Jackson et al 



(|2010l ) found that beach renourishment programs in estuaries can enhance shore protection, 



but can decrease habitat suitability by changing the beach shape creat ing higher berms and 



Jackson et al 



(I2OIOI) 



wider backshores than would occur under natural conditions. However 
did not focus on any species in particular. The controversy about renourishment vs. species- 
abundanc e and other species patterns, has also involved other shoreline dependent taxa. For 



example, (jBrock et al 



20071 ) studied the influence of beach replenishment on sea-turtles find- 



ing a clear d ecline i n the 



restoration. 



nest ing success an d abundance in the sea s on after the anthropic 



Mennl (|2002allbl ); 



Greene! ( 2002 ) 



Guilfovle et al 



(|20061 ): 



de la Huz and Lastra 



(|2008l ) are in agreement that most of the actual beach renourishment projects worldwide 



disturb the food-web structu re of the coas tal hab itat ecosystem impact i ng al l the species 



occupying the affected niche. 



Mennl (|2002allbD . and 



de la Huz and Lastral (j2008l ) noticed the 



evident linkage between the geomorphodynamic structure of the beach and the quality of 
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the habitat in sustaining species (from microorganisms to bi rds). Dredgi ng for beach renour 



ishments was found also to impact coral reef communities (jJaan 



20001 ) . The three factors 



of beach intervention c osts, bi o diver sity protection, and potential income, make the a-priori 



adaptive management (j Thorn 



20001) . ris k assessment, unc ertainty and decision analysis of 



renourishment fundamentally important (jNordstrom 



2005). 



We described in relation to beach renourishment, the role of prey availability, preda- 
tor activity, human activity, and physical features on habitat use and possible survival of 
the Snowy Plover (Charadrius alexandrines) in Florida. For Florida, many technical re- 



(Gore and Chase 




1989 


• 


La 


monte and Douglass 




2002 




Himes et al. 




200G 




US-FWS 




2007; 


Burnev . 


200ol 


USGS-FWS 




2009; 


Lott 




2oogl 


Prunei 




20101). However there 


is still a 


lack of 



quantitative studies addressing the ecological effects of renourishment on shoreline-dependent 
birds. Although the relationship between the use of coastal habitat and s h orebirds has been 
assess e d in much detail by many pr e vious studies, e.g. 



(|2007l ): lHood and Dinsmor 4 |2007j ) 



Sirami et al 



Taft and Haig 



(200S!): lTian 



et al 



(200G); 



(200 



Hood and Dinsmore 



Gan et al 



Colwell et al 



200 



20051 ). the litera- 



as well as for riverine ecosystems in proximity of the coast 
ture about habitat suitability modeling for shorebirds is not extensive. Recently, the Snowy 
Plover has been studied intensively as part of a joint effort of the US Department of Defense, 
the Environmental Protection Agency, and the Dep artment of Energy, for its conservation as 



182 a function o f clima te change and military ac tivity (IChu-Agor et al 



2011 



2010a 



d). 



Converting et al 



2010; 



Converting et al 



(|2010ct ) provided evidence of the fundamental niche of 
SP through a maximum entropy approach, that is composed by the estuarine and ocean 
beaches constituted of alkaline medium/fine white sand and silt 



Convertino et al 



(l2010dl) 



described the source of unc ertainty in data and speci es distribution models for the particular 



case of the SP in Florida. 



Convertino et al 



( 2010a 



found an interesting interannual posi- 



tive feedback between the tropical cyclones in the year prior to a breeding season and the 



SP abundance and range. 



Chu-Agor et al 



Rosa Island along the Florida Panhandle, and 



(2010) modeled the habitat evolution of Santa 



Convertino et al 



(|201lf ) identified a decline in 



the power-law distribution of the habitat patch-size (i.e. the probability of finding suitable 
breeding/wintering habitat patches larger than a given size) for the SP, PP, and RK, through 
coupled modeling of the land-cover and of the habitat suitability as a f unction of the IPCC 
A1B sea-level rise scenario rescaled to 2 m of sea-lever rise. Recently 



Seavev et al 



(120101) 



studied the threat of the Piping Plover as a function of the sea-level rise due to climate 
change in their New York barrier islands habitat. These studies confirm the importance of 
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197 identifying effective interventions, such as renourishments, that support the wildlife needs in 

198 the face of climate change. The purposes of this paper are to: 



199 • Provide a comprehensive overview of the biology of the Snowy Plover in Florida in 

200 relation to other SP populations in the USA, and with other shorebirds in Florida, that 

201 is potentially useful in metapopulation modeling; 

202 • Describe probabilistically the effect of past replenishments on the SP breeding/ wintering 

203 range and abundance, together with the determination of the local-scale geomorpholog- 

204 ical features of the habitat that increase the probability of site-selection and survival. 

205 A comparison is performed with other TER-s species in Florida (PP and RK); 

206 • Assess some management scenarios (specifically "ecologically sustainable" costal restora- 

207 tions) as a function of the environmental cues of the SP to reduce the risk of nest-failure. 
2oa This lays the foundation for subsequent multi-criteria decision analysis for the conser- 
209 vation of the SP. 



2 Materials and Methods 



2.1 Models 



We employed two models: (1) a maximum entropy ap proach model (MaxEnt) (jPhillips et al 



2006 



Phillips and Miroslav 



2008; 



Elith et al 



2010) to quantify the similarities in the habi- 



tat use of the Snowy Plover during the breeding and wintering seasons, and to evaluate the 
habitat preferences of other TER shorebird species; and, (2), a Bayesian approach model to 
evaluate the probability of coarse-scale site-selection for different species after replenishment 
events, sampling from the posterior probability and based on historical data. A detailed de- 
scription of the biological data used in the models is contained in the Supplementary Material 
of the manuscript. 



The ma ximum entropy model adopt ed (MaxEnt) is fully described in ( Phillips et al 



20061 ) and ((Phillips and Miroslav 



20081) . MaxEnt has alre ady been applied i n mod eling 



223 the habitat suitability of the SP during the breeding season by 



Convert ino et al 



Convertino et al 



( 2010d ) and 



(|2010d| ). The estuarine and ocean beaches composed of alkaline medium/fine 

225 white quartz sand and silt were found to be the most suitable for SP during the breeding sea- 

226 son. Here we used the same approach, considering as explanatory variables at the regional 

227 scale (the entire Gulf coast of Florida ) the land-cover and the geology map. The land- 



cover from NO A A (Klemas et al 



19931 ) has been translated into SLAMM land-cover classes. 



SLAMM (Sea Level Affecting Marshes Model) (IClough , 



late the effect of sca-lcvcl rise on the coastal habitat 



20111) . Tables 1 and 2 in 



Convertino et al 



2006) is the model we used to simu- 



Chu-Agor et al 



2010: 



Convertino et al 



(|2010q ) report in detail the land-cover and the 



geology cl asses. The fundamental niche for Piping Plovers and Red Knot is d etermined by 



MaxEnt (jPhillips et al 



2006 



Phillips and Miroslav 



2008: 



Elith et al 



20101 ) for the same 



234 geographical domain of the SP and with the same environmental variables (Figured]). Occur- 

235 rences for the SP, PP, and RK are for the 2006 wintering distribution. The habitat suitability 

236 maps are the average over 30 replicates performed for at least 10,000 random background 

237 points. 



239 To quantify the apparent relationship between beach nourishment projects and SP, PP, 

240 and RK depicted in Figure [3] we used a Monte Carlo procedure to sample from the posterior 

241 distribution of the binomial probabilities that a region is a nesting or wintering ground 

242 conditional on whether or not the same year or the previous year the region experienced a 

243 renourishment event. Table [T] lists the number of nests by year in the breeding and in the 

244 wintering season, the number of renourishment events, and the average number of fledglings. 

245 For example, the data in the gray background in Tablc[T]are those considered for the Bayesian 

246 inference in the seasons 2005-2006 (the renourishment interventions in 2004 are considered 

247 when analyzing the 2005 breeding season). In the years 2008-2010 there were renourishment 

248 events in the three regions considered (Pensacola/Eglin, Tyndall, and Peninsula, see Figure 

249 [I}. Data about the 2008-2010 renourishmcnts arc reported in the Supplementary Material. 

250 The median number of nests per year for the SP, per region is 57. Here we define a nesting 

251 ground as an area having at least 10 nesting sites. A wintering ground is an area in which 

252 at least 2 adult individuals were observed. The threshold of 10 nesting sites is found to 



(Convertino et al. 




2010a) 


2010c 


) the occurrence of 



2010al) . Because the solitary behavior of shorebirds (|Convertino et al 



256 The regions are considered indepen dent because of the fidelity of SP and because of the 



257 small dispersal range of the species (jColwell et al 



Convertino et al 



2007b 



Stenzel et al 



2007 



FWC 



2010; 



2010al) . For PP and RK the three sites can be considered independent 

259 because these shorebirds use the areas as wintering and stopover areas for a limited period 

260 during which the inter-site movement was observed to be very limited. We considered all 

261 the renourishment projects that happened in every region in the previous year or before 

262 the nesting season in the same year, regardless of the size of the renourishment project. 



Specifically we considered all the areas that in the period 2002-2010 were subjected to at least 
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264 one renourishment event. Those areas were considered as the potential breeding/wintering 

265 regions in the Bayesian inference. The occurrence of a nesting or a wintering ground is then 

266 checked in these regions. Let Y be a random variable having a value of one if the region is 

267 not a nesting or wintering ground and zero otherwise and X be a random variable having 
26a a value of one if the region was affected by a beach nourishment in the previous or in the 

269 same year and zero otherwise. Then the odds ratio (OR) of a region not being a nesting or 

270 wintering ground given at least one renourishment event the year before or the same year 

271 relative to the region not being a nesting ground following a year without a renourishment 

272 is given by: 

np = P{Y = l\X = \)P(Y = Q\X = 0) tti(1 - ttq) 

P{Y = 0\X = 1)P(Y = 1\X = 0) 7r (l-7ri)' [> 

273 The likelihood function given by the product of binomial distributions Y± ~ binomial(7Ti , m) 

274 and Yq ~ binomial(7To, no)- Assuming beta priors for the probabilities 7To,7ri with parame- 

275 ters (do, bo) and (ai, b\), our posteriors are given by 7i"o|y ~ beta(yo + no + bo — yo) and 

276 tti \y ~ beta(j/i +oi, n\ + b\ —y\). Assuming a uniform prior on the distributional parameters, 

277 we simulate posterior probabilities directly from the posterior distributions and compute the 

278 odds ratio as the fraction of the odds of an "empty ground" in the breeding season following 

279 a year with at least one beach nourishment project, to the odds of an "empty ground" in the 

280 season following a year without a beach nourishment. The Bayesian approach was performed 

281 for SP in the breeding and wintering season, and for PP and RK for the wintering season. 

282 3 Results and Discussion 

283 As for the TER-s shorcbirds analyzed, the Tyndall area hosted about 40 %, 25 %, and 22 

284 % of the SP, PP, and RK populations. Additionally for the Wilson's Plover (competitor of 

285 the SP) 83 % of the population was in the Peninsula in the breeding season, confirming the 

286 importance of those Florida west coast Gulf beaches. The high presence of the WP and SP in 

287 those areas can explain the high favorability of the Panhandle for these shorcbirds during the 

288 breeding season. Pensacola and Eglin areas hosted on average 8-10 % of the SP population 

289 in the winter and in the summer season. For PP and RK those areas represent less than 1 

290 % of the population. The Peninsula and the Atlantic coasts are the main wintering grounds 

291 for the migratory PP and RK. The PP Peninsula population was 38 %, and the Atlantic 

292 population was 33 %. For the RK the proportions were 55 % and 20 % for the Peninsula and 

293 the Atlantic. Broad-scale estimates of the hedge-rate for Snowy Plovers nesting in Florida 

294 was estimated from data (Table [!}, however juvenile survival rates and adult survival remain 
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295 still unknown. According to the trend of the average number of fledglings, an increase in 

296 the breeding population of SP is expected, supposedly because the increasing care in their 

297 conservation. 



Figure [5] reports the suitability index maps (SI) for the wintering season of SP (Fig. [3J 
a), PP (Fig. [5J b), and RK (Fig. [51 c) in the Panhandle-Big Bend- Peninsula region. We 
decided to model PP and RK in the same geographic domain where the range of the SP oc- 
curs in order to perform a comparative analysis for the habitat use of the studied shorebirds. 
The constraints of the habitat suitability model (MaxEnt) are the adult-pairs occurrences 



tat classes ( 


Chu-Aecor et al. 




2010 


■ 


Convertino et al 


(Convertino et al. 


2010c 


), at a resolution of 120 m. 



2010d ). and the geology GEO classes 



calculated the probability map (fr om to 1) assum i ng a r egularization parameter equal to 



one, pseudoabsences placed as in (jConvertino et al 



2010cJ), and 25% of the occurrences as 



trainin g sample. The breeding habitat suitability map for SP is reported in (jConvertino et al 



2010c). From the habitat suitability map that can be considered the predicted fundamental 
species distribution, we attributed the class Suitability Index SI=60 to every pixel whose 
probability values are > 0.2 (threshold value) of the species distribution model. In this 
way considering only the pixels f or SI > 60 the g e ograp hic range for the shorebirds ana- 



lyzed is successfully reproduced (|Convertino et al 



20111) . All the pixels for SI < 60 are 



categorized as unsuitable. SI=60 is considered the lowest score associated with consistent 
use and breeding/ wintering, SI=80 is the score typically associated with successful breed- 
ing/persistent wintering, and SI=100 is f or the best habitat with the highest surviv al and 



reproductive success or stable wintering (jMaika et al 



2007 



Convertino et al 



20111) . Fig- 



ure^ (d, e) shows the response function or conditional probability of presence as a function 
of the two explanatory variables, the land cover and the geology. The logistic prediction 
changes as each environmental variable is varied, keeping all other environmental variables 
at their average sample value. In other words, the response curves show the marginal effect 
of changing exactly one variable, whereas the model may take advantage of sets of variables 
changing together. From the response curves it is possible to note the similar habitat prefer- 
ences of SP in the breeding (dashed blue line) and wintering (dashed red line) seasons. In the 
winter the SP se em to use more the ocean bea ches (class 12) than in the breeding season as 



documented by (jLamonte and Douglass 



20021). The Piping Plover has a habitat preference 



similar to that of the Snowy Plover. Results of Figure [2] (b) (green curve) show that the PP 
also occupies scrub/shrub transitional marsh and salt marsh areas with higher probability 
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than the SP (class 7, 8, and 9 respectively) as confirmed by observations 



Elias et al 



(2000) 



Elliott Smith and Haia (|2004h . Results also show that the PP seems to use more the ocean 
beach than SP does in the winter (higher P(X — SLAMM=12)). The RK in the winter seems 
to prefer more estuarine beaches (class 10). However, relatively high values of the probability 
of occurrence are observed also for the other SLAMM classes in comparison with SP and 
PP. The presence of medium/fine alkali sand and silt is less a requirement for the RK. RK 
is very adaptable to any substrate and in particul ar the suitability is high also for peaty- 



substrate habitats (class 12 GEO) as reported by (|N 



les et al 



2008). The minimization of 



Convertino et al 



2010dh . for example 



33s the uncertainty of these results has been obtained in 

339 considering the positioning error of recorded occurrences and the spatio-temporal gaps be- 

340 tween occurrences and land-cover maps. The differences in the habitat use among shorcbirds 

341 in the winter underline the importance of careful restoration planning policies that try to 

342 accommodate the needs of all the sensitive species. From the habitat suitability modeling, 

343 it appears that the resident SP is the most sensitive of the three species in relation to the 

344 habitat use (ocean and estuarine beaches) when subjected to variations due to renourshimcnt 

345 events. Piping Plovers have a very similar habitat use in the wintering season when they 

346 migrate to the coast of Florida; however they seem to be more resilient to habitat variations 

347 than Snowy Plovers because of their wider habitat preferences. Red Knots are the least sen- 

348 sitive species in relation to renourishment projects that modify the estuarinc/ocean beaches, 

349 since they show the broadest spectrum of habitat preferences. 



Figure [3] reports the observed distribution of SP nests and SP, PP, and RK adult-pairs 
by year for the Pensacola/Eglin, Tyndall, and Peninsula study areas. The renourishments 
in the 2005-2006 are represented in Figure however in the analysis the whole 2002-2010 
period of data is considered (Table [1]). In each plot the breeding and wintering distributions 
of SP, PP, and RK are reported in the same year or in the year following the replenish- 
ment events. Only few nesting and wintering grounds occur in the locations where beach 
renourishment events occurred the previous year or the same year. In the data available 
there is not a complete information about the detailed timing of each intervention. As a 
result it is not possible to fully understand if the renourishments were performed during the 
wintering season or, less likely, during the breeding season and what was the duration of 
the project. However this information was partially compiled using the renourishment data 
reported in the Supplementary Material. The purpose of the current study was to infer the 
feedback between renourishment and TER-s shorebirds considering their spatial occurrences. 
Consequently we can not distinguish the direct or long-term effect of the renourishment on 
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365 shorebird species. Renourishments events occured also during the beginning of the SP breed- 

366 ing season along the Florida Panhandle. Those interventions surely disturbed directly the 

367 site-selection processes. SP tend to have high-fidelity to nesting sites. Moreover, they are soli- 

368 tary nesters, unlike colonial watcrbirds, and direct disturbances during the breeding season 

369 can seriously increase the fragmentation of the suitable habitat that affects their distribution 

370 and nesting success. The integrity and extension of a species' habitat are features that af- 

371 feet the survivability of the single individuals and the extinction risk of the whole population. 

372 

373 Figure [4] show the posterior probability of absence P(A > a) of the odds ratio of not 

374 being a SP breeding ground, or a SP, PP, and RK wintering ground, based on the historical 

375 data. For SP over the years 2002-2010 we have one case (2004) in the Panhandle in which 

376 there was not any renourishment activity and there were more than two nesting grounds. In 

377 the same year (2004) for the Peninsula there were 2 renourishment projects and no nesting 

378 grounds. A value of P(A > a) above one indicates a relatively higher probability of a region 

379 not being a nesting ground following a year with a renourishment. The distribution is skewed 

380 to the right with a mode of about 1.7. The median value of the odds ratio indicates that it 

381 is 2.5 times more likely that a region will not be a nesting ground for a SP following a year 

382 with a replenishment event (dashed solid curve). A 90% confidence interval for the odds ratio 

383 is (2, 30). Specifically, for the SP we have y\ = 12 cases of SP nesting grounds (as regions 

384 having at least 10 nest counts) over a sample of ri\ = 46 breeding seasons (counting separately 

385 all the renourished areas in 2002-2010) that were exposed to a beach nourishment the year 

386 before or the same year, while we have yo = 30 cases of nesting grounds over a sample of no 

387 = 45 breeding seasons that were not exposed to a previous year/same year renourishment. 

388 Assuming a uniform prior (ao = ai = £>o = b\ = 1), we simulate 10 4 Monte Carlo posterior 

389 probabilities directly from the posterior distributions f{^o\y) and f{^i\y) and compute the 

390 OR as given in Eq. [1] Statistics on the OR are derived directly from the posterior samples. 

391 The results for the historical wintering distribution indicate that it was 2.3, 3.1, and 0.8 times 

392 more likely that a region was not a wintering ground following a year with a replenishment 

393 event for the SP (solid red curve), PP (solid brown curve) and RK (solid green curve) respec- 

394 tivcly. The fit of the calculated probability is made by a lognormal distribution with different 

395 shape parameter. The posterior probability of absence for the breeding season is stable for 

396 values of the cutoff in considering an area a breeding ground of 10 ± 4 nests . The r e hav e 



been speculations that biological mechanisms induce log-normal distributions (jKoch 



1966). 



In the majority of plant and animal communities, the abundance o f species follows a trun 



cated log-normal distribution (jSueihara 



1980 



Limpert et al 



200 lh . Our conclusion is that 
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regardless of the lack of detailed information about renourishment projects we are confident 
in assessing the negative feedback between past renourishment projects and SP, PP, and RK 
at the local-scale and at the macroscale for the whole region considered. The RK appears to 
be the least affected shorebird by beach renourishmcnts. Our goal is to emphasize ecological 
sustainable restoration projects that take into account the habitat preferences of resident and 
migratory species. For example, the use of submerged geotexile groins that trap the sand 
ncarshore, has no impact on the existing beach/ocean habitat communities and increases 
the beach area naturally. A successful application of geotexile groins has been tested in the 
Charlot te county shores in Flori da also during the devastating 2004-2005 tropical cyclone 



season (Beach Restoration 



2005) 



4u The large-scale prediction of the habitat suitabi lity for the whole SP Florida range 



412 (jChu-Aeor et al 



2010: 



Convertino et al 



2010c . 



20111) coupled with the local niche analysis 



as a function of physiognomic and biotic variables (here ass umed to be invariant to c limatic 



414 changes), environmental variables (e.g. tropical cyclones (jConvertino et al 



2010ah ). and 



human interventions (e.g. renourishments) will significatively help the estimation of the ex- 
tinction and decline risk of Snowy Plovers in population viability models. These quantitivc 
studies will enhance the adoption of effective conservation policies, e.g. through multi-criteria 
decision analysis, for the conservation of imperiled species. 



■119 4 Conclusions 

420 The identification of habitat cues is critical for appropriate ecosystem management aimed to 

421 the conservation of rare and declining species. A comprehensive Hutchinsonian description 

422 of species is a necessary step in accomplishing this goal, combined with the understanding of 

423 the species's biogeographical distribution. Local scale analysis and macroscale studies need 

424 to be combined to adopt efficient conservation policies. Here we focused on the relation- 

425 ship between resident and migratory Threatened, Endangered and at-Risk (TER) shorebirds 

426 and renourishmcnts along the Florida Gulf coast. In particular the focus was on the Snowy 

427 Plover, a resident state-designated threatened shorebird. In comparison we analyzed the Pip- 

428 ing Plover, a migratory federally-designated threatened shorebird and the Red Knot which 

429 is a least-concern shorebird in Florida. The analysis was performed on nest occurrences and 

430 not on nesting success or chick survival data, which are not available. However the observed 

431 spatio-temporal correlation between renourishmcnts and nest/adult-pairs of shorebirds can 

432 provide insights into the causes of the reduced number of nests after renourishment interven- 
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433 tions. The following conclusions are worth mentioning. 



434 • Based on the 2006 wintering and breeding counts Tyndall is confirmed to be the hotspot 

435 of the shorebird species richness in Florida both for resident and migratory birds. The 

436 differences in the breeding and wintering distribution of the SP population (in the Pan- 

437 handle 80 and 60 % in the breeding season and in the wintering season; in the Peninsula 

438 20 and 40 % respectively) and the almost unaltered total pairs count can potentially 

439 confirm the predicted movement of SP (~ 20%) to the lower and warmer latitudes of 

440 the Peninsula beaches during the winter. However SP from other Gulf state beaches 

441 have been sporadically observed in Florida. The result confirm the observations fo r 

442 the Snowy Plover subpopulation in the West coast of the USA (jStenzel et all 119941) , 

443 in which the inter-seasonal dispersal happens for longer distances than within seasons. 

444 Renourishment interventions need to carefully preserve Tyndall as focal area for shore- 

445 birds richness. Moreover dispersal patterns have to be considered in order to reduce 

446 the potential direct disturbance of renourishmcnts. Nonetheless further studies about 

447 the shorcbirds movements patterns in Florida are needed; 

448 • During both nesting and brood-rearing stages of breeding, Snowy Plover selection of 

449 habitat and productivity are influenced by a combination of abiotic and biotic factors 

450 including human disturbance, predator abundance, prey availability and the physical 

451 features of the habitat. However, we did not study the factors influencing different 

452 breeding stages due to the lack of data. The habitat su i tability model based on the 

453 princ iple of maximum entropy (MaxEnt) (jPhillips et all l2006t iPhilhps and Miroslav . 

454 120081 )) did not observe consistent differences between the habitat preference of SP in 

455 the breeding and in the wintering season. In the winter the SP seems to utilize more the 

456 ocean beaches than in the summer, possibly because the brood-rearing habitat is less 

457 used than during the breeding. In the wintering season the habitat use of the migratory 

458 PP is very similar to that of the SP. In contrast, the RK seems to have a larger spectrum 

459 of habitat preferences. The PP utilizes the ocean beaches more than the SP. However, it 

460 also prefers scrub/shrub transitional marsh and salt marsh areas. The RK is observed 

461 occasionally also on peaty-substrate banks. We did not consider any possible direct 

462 interaction (such as interspecies density-dependence) between the shorebird species. 

463 PP and RK seem more resilient to the effects of renourishment projects (that modify 

464 estuarine and ocean beaches) that are not focused on the conservation of the wildlife 

465 habitat because their habitat preference is wider than SP; 

466 • We argue that the decrease in the documented nest abundance of SP resulted from 
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an altered cross-sectional beach profile which is not favorable for nesting and foraging. 
As a consequence the nesting success is reduced because of habitat modifications and 
possible direct disturbance of renourishment interventions. The profile subsequently 
improved in later seasons as the beach equilibrated to a more natural slope and surface 
roughness more significantly. A negative feedback was found between SP, PP, and RK 
and historical renourishment projects. Results based on a spatio-temporal Bayesian 
inference show that it was 2.1, 3, and 1 times more likely that a region was not a 
wintering ground following a year with a replenishment event for the SP, PP and RK, 
respectively. Despite the fact that the inference for the TER shorebirds in the winter 
is based only on the 2002 and 2006 census we believe it is significant and in agreement 
with the observed behavior. Considering the census of nine breeding seasons and the 
renourishment events in the same period, the median value of the odds ratio indicates 
that it was 2.5 times more likely that a region was not a nesting ground for a SP in 
the same or in the season following a renourishment event. The higher median of the 
OR in the breeding season indicates that the renourished beaches were more altered 
in the summer season immediately after the renourishments (performed mostly in the 
winter). The physiochemical properties of the dredg ed material for the renour i shment 



were e quivalent to those of the existing in-situ sand (jLott 



2009: 



US ACE . 



2009: 



BBC'S 



2010bl ). It can therefore be concluded that the ecogcomorphological alteration of the 
habitat was the only cause of the diminished occurrence of shorebirds in the replenished 
sites. Natural processes such as ovcrwash and tropical cyclones shaped the renourished 
beaches, bringing them to their non-altered configuration, for e.g. with the creation of 
ephemeral pools and bay tidal flats that constitute the favorite habitat for SP and PP 



( Converting et al 



2010ah : 



• We emphasize the importance of an a-priori planned "ecological sustainable renourish- 
ment" (e.g., performed by submerged gcotexile groins that preserve the beach cross- 
section profile) that consider the triality among dredging costs, protection of the coastal 
structures, and the potential income from the value of the preserved biodiversity and 
the enhanced recreational activities on the extended beaches. Renourishment projects 
designed to create high quality brood-rearing habitat characterized by sparse low veg- 
etation, and the maintaining of high-prey foraging habitat is an important part of 
shorebirds conservation. 
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Table Captions 



Table 1. Breeding (6) (iFWd l2010l: lAlliancd . 120101) and wintering (w) (|USGS-FWSl . 



20091 ) counts (SP nest and adults pairs respectively) for the Panhandle, Peninsula and the 



whole populations in the seasons 2002-2010. A SP pair is con sidered to exist for every nest 



20101) . r (P and T stand for Pen- 



744 count. The number of renourishment events (jBSRC . 

745 sacola/Eglin, and Tyndall), refers to the renourishments made the same year or the year pre- 

746 vious to the breeding and wintering season. We considered as potential breeding/ wintering 

747 regions in the Bayesian inference all the areas that in the period 2002-2010 were subjected to 
74a at least one renourishment event. The number of fledglings is estimated from the observed 

749 counts. The data in the gray background arc those considered for the Bayesian inference in 

750 the seasons 2005-2006 (Figure |3J). 



Figure Captions 

Figure 1. (a) Map of the Panhandle-Big Bend-Peninsula (PBBP) study area along the 
Gulf coast of Florida, and closeups of the three focal study sites (Pensacola, Eglin, and Tyn- 
dall (Apalachee Bay)) in which there is a high density of military areas, federal reserves and 
state parks. The red dots re fer to t he 2006 Snowy Plove r (SP) breeding census performed 



757 along the Florida coastline (|FWC 



2010: 



SP adult-pairs in the winter season (USGS-FWSl. 



Wilson's Plovers (FWC 



2010: 



Alliance 



Alliance 



I 201oi 



20101) . The gold dots are the observed 



2009) . The blue dots are the nests of the 



(see Supplementary Ma terial). The c ritical 



2010al) . (a), 



beaches are depicted in red based on the 2009 Critical Erosion Report (|BBCS , 
(b), and (c) are the 2010 satellite images of the selected study sites with the US military 
bases involved in the study deline ated in red, (b) w intering distribution of Piping Plover 



(PP) and Red Knot (RK) in 2006 l|USGS-FW£ , 



2009) 



Figure 2. Suitability Index maps derived from the average over 30 habitat suitability 
realizations calculated by MaxEnt for the Snowy Plover (SP) (a) , Piping Plover (PP), an d 



2010c). 



Red Knot (RK) as a function of the land-cover, the geology layers (| Converting et al 
and the 2006 winter adult-pairs occurrences. The conditional probability P{X\Y) to find a 
nest or an adult-pair is plotted as a function of the continuous explanatory v ariable Y at reso 



lution 120 m, that is the land-cover translated into SLAMM habitat classes (|Chu-Agor et al 



2010 



Convertino et al 



2010a ) (d), and the geology GEO (e), for the model run keeping all 
other environmental variables at their average sample value. X is a SP nest in the SP breed- 
ing region or a SP, PP, and RK adult-pair in the wintering season. 



775 Figure 3. Distribution of SP nest sites (dots) and SP, PP, RK adult-pairs sites (circles 

776 proportional to the adult-pairs abundance) by year for Tyndall (a), Pensacola/Eglin (b), and 

777 Peninsula areas (c, d). The renourishments in the 2005-2006 period are represented. In each 

778 plot is reported the breeding and wintering distribution of SP, PP, and RK in the same year 

779 or in the year following the replenishment events. In the b ackground the habitat su itability 



2010d). A dot 



7so is represented for the buffer of 10 km from the coastline ([Converting et al. 

781 may represent more than a single nesting site. Within the plot the text indicates the renour- 

782 ishment R in the previous year that is represented by a continuous light-blue line in the map. 

783 Red arrows indicate the sites with positive feedback SP-renourishment. 

784 

785 Figure 4. Posterior probabilities of absence P(A > a) of the odds ratio for SP in the 

786 breeding season, and SP, PP, and RK in the wintering season. The odds ratio is the ratio 

787 of the odds of a nesting or wintering ground in the spring following a year with at least one 

788 renourishment event to the odds of a nesting or wintering ground in the spring following a 

789 year without a renourishment intervention. For the breeding SP the median odds ratio is 

790 2.5 and the mean is 4.9. For the wintering SP, PP, and the RK the median odds ratio is 

791 2.3, 3.1 and 0.8 respectively. The maximum likelihood estimate is a lognormal distribution 

792 with different values of the shape parameter (histogram only for SP) and the coefficient of 

793 determination for SP, PP, and RK is on average R 2 = 0.92. 



25 



Tabic 1: 











Pairs 












Year 




Panhandle 


Peninsula 




Whole 


Fledglings 




b 


w 


r 


b 


w 


r 


b 


w 




2002 


128 


228 


1 (P) 


65 


103 


1 


193 


332 




2003 


9 




1 (P) 


1 




4 


10 






2004 


57 










2 


57 






2005 


4 




2 (P, T) 


93 




6 


97 




1.082 


2006 


235 


175 


2 (P, T) 


68 


137 


7 


303 


312 


1.075 


2007 


48 




1(P) 


96 




6 


144 




1.524 


2008 


394 




1(P) 


78 




1 


472 




2.000 


2009 


1051 




2(P) 


195 




1 


1246 




0.860 


2010 1 


241 




2 (P, T) 


66 




5 


307 




1.390 



*Data updated to the 30 July 2010 of the breeding season. 
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